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' We investigate the gravitational energy emission of an ultrarelativistic particle radially falling 

^ ' into a D- dimensional black hole. We numerically integrate the equations describing black hole 

gravitational perturbations and obtain energy spectra, total energy and angular distribution of the 
emitted gravitational radiation. The black hole quasinormal modes for scalar, vector, and tensor 
perturbations are computed in the WKB approximation. We discuss our results in the context of 
black hole production at the TeV scale. 
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Brane-world models describe the visible universe as a four-dimensional brane embedded in a higher-dimensional 
bulk 1]. A generic consequence of the brane-world scenario is that the fundamental gravitational scale is lower than 
> , the observed Planck scale. In some models, the fundamental scale is lowered to values that would be accessible to next- 
generation particle colliders, thus enabling laboratory-based studies of strong gravitational physics via perturbative 
[ 13 and nonperturbative events |3| . Ultrahigh energy cosmic rays could also probe trans-Planckian energies i. The 
• possibility that strong gravitational effects such as black hole (BH) and brane formation could be observed in the near 
^ ' future has sparked a lot of interest in the investigation of nonperturbative gravitational phenomena in hard-scattering 
• ^ ■ events [3. (For a review and more references, see Ref. 0). 

■ Trans-Planckian BH formation at energy scales much larger than the fundamental gravitational scale is a classical 
process The event is dominated by the s-channel and the initial state is modelled by two classical shock waves 
with given impact parameter. In this context, a major issue is the estimate of the coUisional energy loss. The hoop 
conjecture states that the collision of two particles ij with center-of-mass (cm.) energy Ecm and impact parameter 
smaller than the Schwarzschild radius ?'s(i?cm) forms a trapped surface f^. This event is formally described by the 
process ij — > BH-|-£'(X), where E{X) denotes the coUisional "junk" energy which does not contribute to the BH mass. 
The junk energy includes a bulk component of gravitational radiation and other possible non-standard model gauge 
fields, and a brane component of standard model coUisional by-products carrying the charge of the initial particles. 
The newly formed BH is expected to decay first by loss of gauge radiation into the bulk and then by thermal Hawking 
emission. The Hawking evaporation ends when the mass of the BH approaches the fundamental gravitational scale. 
At this stage the BH either decays completely by emitting the residual Planckian energy or leaves a stable remnant 
with mass about the Planck mass . Most of the observable signatures of BH formation come from Hawking's phase 
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and strongly depend on the initial BH mass '9|. Hence, a precise calculation of the coUisional energy loss is essential 
to the phenomenology of BH formation. 

A numerical estimate of the total collisional energy loss for spherically symmetric BHs in Z? > 4 dimensions has 
been given by Yoshino and Nambu (YN) lo] (see also Ref. I^)- The YN approach evaluates the total junk energy 
E{X) by investigating the formation of the BH apparent horizon The colliding particles are assumed massless, 
uncharged and pointlike. Each particle is modelled by an infinitely boosted Schwarzschild solution with fixed energy. 
This solution describes a plane-fronted gravitational shock wave corresponding to the Lorentz-contracted longitudinal 
gravitational field (Aichelburg-Sexl wave) p"s |. The collision is simulated by combining two shock waves travelling in 
opposite directions. The apparent horizon arises in the union of the two shock waves. The junk energy is estimated 
by comparing the initial cm. energy to the BH mass. The result is that the collisional energy loss depends on the 
impact parameter and increases as the number of spacetime dimensions increases. 

The YN method allows estimation of the total junk energy in the classical uncharged point-particle approximation. 
However, it cannot discriminate between different components of E{X), which is theoretically and experimentally 
most important. In a realistic BH event such as a proton-proton collision at LHC the BH is formed by the 
collision of two partons. The bulk component of the junk energy is dominated by gravitational radiation and is 
invisible to the detector. The gravitational junk energy and the invisible component of Hawking emission (neutrinos, 
gravitons . . . ) add to the total missing energy of the process. Therefore, the knowledge of collisional energy loss in 
gravitational emission should provide a good estimate of the different sources of energy loss and missing energy. 

An accurate estimate of the gravitational collisional energy loss would require the use of the full non-linear Einstein 
equations in D dimensions. This is a formidable task, even in four dimensions. Recently, significant advances in 
numerical relativity allowed stable numerical simulations of BH-BH collisions for initial BH separation of a few 
Schwarzschild radii in the non-linear Einstein theory. The gravitational waveforms predicted by these simulations 
are in excellent agreement with analytical results from first and second order perturbation theory |l5j |. Since the 
linearization of the Einstein equations yields results which are surprisingly close to the full theory (see, e.g., Ref. |l6j|). 
BH perturbation theory is likely to provide accurate estimates of gravitational wave emission in higher-dimensional 
spacetimes. Relying on this result, we compute the gravitational wave emission in higher dimensions via a perturbative 
approach. Our computation is the first of this kind to our knowledge. 

The formalism for the computation of gravitational wave emission from perturbed BHs was developed by Regge 
and Wheeler ^3 ^^'^ Zerilli Ts'l, who reduced the problem to the solution of two Schrodinger-like equations. Davis 
et al. [l9| computed the energy radiated in the radial infall of a particle of mass mg starting from rest at infinity 
into a four-dimensional BH of mass Mbh ^ ?7io- This study was later generalized to the radial infall of a particle 
with finite initial velocity or starting at finite distance from the BH |20I. (F or a more comprehensive introduction to 
BH perturbation theory see, e.g., Refs. r21]). Cardoso and Lemos |22L |23| have recently investigated the plunge of 
ultrarelativistic test particles into a four-dimensional static BH and along the rotation axis of a Kerr BH, improving 
early estimates by Smarr p3 |. In this paper we generalize these results to higher dimensions by computing the 
gravitational radiation emitted by an ultrarelativistic particle falling into a D-dimensional spherically symmetric BH. 
Since wave propagation in odd-dimensional curved spacetimes is not yet fully understood, we restrict our investigation 
to even dimensions. (Wave late-time behavior and propagation are very different in odd- and even-dimensional 
spacetimes [l^ll^ll^- Moreover, open issues in the definition of asymptotic flatness do not allow an unambiguous 
definition of "gravitational waves radiated at infinity" in odd dimensions.) 

We model the particle collision as a relativistic test particle plunging into a BH with mass Mbh — Ecm- We use 
recent results of D-dimensional gravitational-wave theory by Cardoso et al. and the D-dimensional extension of 
Zerilli's formalism by Kodama and Ishibashi (KI) |29l ISOl I3lll32 | . which reduces the problem to the solution of three 
Schrodinger-like equations. Our method provides a simple and relativistically consistent estimate of the collisional 
gravitational emission in higher dimensions. We derive the emitted energy in terms of the wave amplitude and study 
the angular dependence of the radiation using the KI formalism. We also present a systematic calculation of BH 
quasinormal modes (QNMs) for the different perturbations in the WKB approach, extending recent calculations by 
Konoplya ^33, ^sil • We show that there is a significant relation between the QNM frequencies and the spectral content 
of the emitted radiation. 

The outline of the paper is as follows. In Sec. |21 we introduce our notations and the basic equations. In Sec. El 
we briefly describe our numerical approach to the computation of gravitational wave emission (details are in the 
Appendices). Section ^contains the main results of the paper. Conclusions are presented in Sec.|31 

2. PERTURBATION EQUATIONS AND QUASINORMAL MODES 

In the next subsection we introduce the background metric and the KI perturbation equations 's^. In Subsec. 12.21 
we describe the method to compute the BH QNMs. 
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2.1. Background metric and perturbation equations 

The spherically symmetric BH in D — n + 2 dimensions is described by the Schwarzschild-Tangherlini metric ^35.J 

ds' = -fir)dt' + ^+r^dnl, (1) 
f(r) 

where dQn is the metric of the n-dimensional unit sphere iS", and 

The BH mass A/bh is given in terms of the parameter M by 

_ nMAn 
JWbh - z — ^ , 

where An = 27r("+i)/Vr[(n + l)/2] is the area of S'", G„+2 is the (n + 2)-dimensional Newton constant, and c is the 
speed of light. We will set Gn+2 = 1 and c = 1 in the following. The {n + 2)-dimensional tortoise coordinate is 
defined by 

' (4) 



Integrating Eq. we find 



dr f{r) 



2M ln{r/a, - 1) 



where 



a, = (2Af)i/("-i)e2"'^/("-i\ (j = 0, . . . , ti - 2) , (6) 



and the integration constant has been chosen to make the argument of the logarithm dimensionless. Here and 
throughout the paper we use the notations of Refs. 29, 30]; the indices (/i, v), and (a, b) denote the coordinates 

of the D-dimensional spacetime, the coordinates of 5", and the coordinates of the two-dimensional spacetime (t, r), 
respectively. 

Kodama and Ishibashi |30l | showed that the gravitational perturbation equations for this metric can be reduced to 
Schrodinger-like wave equations: 

'^^ uj^-v]^ = 0, (7) 



where the potential V depends on the kind of perturbation. Setting x = 2M/r"^^, the potential for scalar perturba- 
tions is 

y f{r)Q{r) 

- -j^^ , (8) 



where 



and 



n(n -\- 1] 

= l(l + n-l), 1 = 0,1,2,... m = K^-n, H{r)=m+-— — '-x , (9) 



Q(r) = n^{n+lfx^+n{n + l)[A{2n'^--in + 'i)m + n{n-2){n-A){n + l)\x^ 

- I2n[{n - A)m + n{n + l){n - 2)]mx + IQm^ + An{n + 2)rr? . (10) 

Equation (jSJ reduces to the Zerilli equation T^l for n = 2. The potential for vector perturbations is 



/ 



n 



{n - 2) 3n^M 



2r"' 



(11) 
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where 

kI, = l{l + n-l)-l, 1^1,2,.... (12) 
Equation l(TT]l reduces to the Regge- Wheeler equation ^ for n = 2. FinaUy, the potential for tensor perturbations 



Vt = ^ 



n(n~2) n^M 

— H 

4 2r"-i 



where 



k| = /(? + n - 1) - 2, 1 = 1, 2, . 



(13) 



(14) 



Equation (|13|l was derived by Gibbons and Hartnoll |3g in a more general case (see also [s^, where a Gauss-Bonnet 
term is included) and has no equivalent in four dimensions. 



2.2. Quasinormal modes 



The knowledge of the QNM frequencies of multidimensional BHs enables a clear physical interpretation of their 
gravitational emission. QNMs are free damped BH oscillations which are characterized by pure ingoing radiation at 
the BH horizon and pure outgoing radiation at infinity. The no-hair theorem implies that QNM frequencies depend 
only on the BH mass, charge and angular momentum. Numerical simulations of BH collapse and BH-BH collision 
show that, after a transient phase depending on the details of the process, the newly formed BH has a ringdown phase, 
i.e., it undergoes damped oscillations that can be described as a superposition of slowly damped QNMs (modes with 
small imaginary part). Furthermore, the QNMs determine the late-time evolution of perturbation fields in the BH 
exterior (for comprehensive reviews on QNMs see Refs. [3^1. 

Gravitational radiation from four-dimensional astrophysical BHs is dominated by slowly damped modes. In the 
following we show that these also dominate the emission of gravitational radiation in higher dimensions and determine 
important properties of the energy spectra. Recently, Konoplya computed slowly damped QNMs of higher-dimensional 
BHs js^, 13 using the WKB method. This method is known to be inaccurate for large imaginary parts, but it is 
accurate enough for the slowly damped modes which are relevant in our context. Therefore, QNM frequencies for 
scalar, vector, and tensor gravitational perturbations are computed here in the WKB approximation. Our results are 
in good agreement with those presented by Konoplya in Ref. 34] (modulo a different normalization). At variance 
with Ref. '34'!, we concentrated on uncharged black holes in asymptotically flat, even dimensional spacetimes. We 
extended Konoplya's calculation in two ways: i) in addition to the fundamental QNM we also computed the first two 
overtones; ii) we carried out our calculations for a much larger range of values of I (Ref. [SJ] only shows results for 
l = 2&ndl = 3). 

The method consists in applying the WKB approximation to the potential in Eq. |(7J with appropriate boundary 
conditions. The result is a pair of connection formulae which relate the amplitudes of the waves on either side of the 
potential barrier, and ultimately yield an analytical formula for the QNM frequencies (for details see Refs. |3^E3|)- 
The WKB QNM frequencies uP are given in terms of the potential maximum Vq and of the potential derivatives at 
the maximum by 



2Vi' (1 



J =0,1,2, 



where 



A = 




lOOa^ 



(15) 



(16) 



(17) 



V" 
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a = J + 1/2 and j is the mode index. The QNM frequencies for the scalar, vector, and tensor potentials of Sec.[5]and 
various dimensions are shown in Tables [TITVI and will be discussed in Sec.^ Let us stress that the application of the 
WKB technique is questionable in a few higher-dimensional cases; for / = 2 and I = 3 the vector and scalar potentials 
in £) > 6 are not positive definite and/or display a second, small scattering peak close to the BH horizon. An accurate 
analysis of these potentials would require a refinement of the standard WKB technique, which is not presented here. 
These special cases are denoted with italic numbers in Tables ITlIIIVI 

We mention that highly damped QNMs of four- and higher-dimensional BHs have recently become a subject of great 
interest in a different context. A few years ago. Hod prmxised to use Bohr's correspondence principle to determine 
the BH area quantum from highly damped BH QNMs |4J. Hod's proposal is quite general: the asymptotic QNM 
frequency for scalar perturbations of a non-rotating BH in t) dimensions is the same as in four dimensions • Quite 
notably, this result holds also for scalar, vector and tensor gravitational perturbations [i^ |43| . Ref. contains a 
partial list of references on recent developments in this field. 



3. INTEGRATION METHOD 



The computation of the gravitational wave emission of an ultrarelativistic particle plunging into a BH requires the 
numerical integration of the inhomogeneous wave equation for scalar gravitational perturbations. (Vector and tensor 
gravitational perturbations are not excited by a particle in radial infall.) The source term 5*'"' for the corresponding 
wave equation in n + 2 dimensions can be calculated from the stress-energy tensor of the infalling particle. Details of 
the derivation are in Appendix 1X1 

The integration in (n + 2) dimensions proceeds as in four dimensions l2r)| |. A good summary of the integration 
procedure can be found in Ref. In this section we simply stress the differences between the four- and the 

(n -|- 2)-dimensional cases. For the sake of simplicity, in our numerical integrations we set the horizon radius = 
(2M)^/("^^' = 1. The equation for the scalar perturbations is 

{£^^+u--Vs)^ = S^-K (18) 

The general solution of Eq. H18() is obtained via a Green function technique as follows. Consider two independent 
{left and right) solutions of the homogeneous equation with boundary conditions <&l ^ e""^''* for — > — oo, and 
~ e"^''* for r* +oo. For r, +oo the left solution is a superposition of ingoing and outgoing waves of the form 



$L - B{uj)e''^''' + C{Lu)e- 



(19) 



The Wronskian is given hy W — 2iLuC{Lu). The wave amplitude is obtained from a convolution of the left solution 
with the source term 



1 



+ 00 



(20) 



The energy spectrum can be expressed in terms of the wave amplitude as (details of the derivation are given in 
Appendix IB|i 



dE 



(21) 



duj 16n n 

where = l[l + n — 1). The Wronskian for a given value of u; is obtained by integrating the homogeneous equation 
from a point located as close as possible to the horizon, and expanding $l as 



where 



O'n+2 



$L~e-'"''- [l + a„+2(r-l) + ...] , 
-(/4 + 2/3 - ;2 _ 2/ + 3) + (n - 2) [-2/3 + P + [n^ + 1)1 - (n^ + An^ + n + 6) /4] 



(22) 



(23) 



(2iw - l)(/2 + / + 1) + (n - 2) [-P + {2iuj - n)l + (n + l)iuj - [n^ + l)/2] 
C{ijj) (and W) can be obtained with good accuracy by matching the numerically integrated to the asymptotic 



expansion 



$7 



a„+2(t^) h.n+2{'^) 



^ ^ a„+2(-t^) ^ 6n+2(-t^) ^ 



(24) 
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where the leading-order coefficient is 



\[P + {n-l)l + n{n-2)/A] 
a„+2(w) — . (25) 



For given n, I and w, the error on the Wronskian and on the energy spectrum is typically of the order of 0(10 



4. RESULTS 



The main results of our work are the computation of the QNM frequencies in the WKB approximation, the 
computation of the energy spectra, and the estimate of the total energy and angular distribution of the radiation 
emitted during the plunge. These results are discussed in detail below. 



4.1. Quasinormal frequencies 

The WKB QNM frequencies for different even values of n are listed in Tables IIIIVI Each line shows the first 
three quasinormal frequencies {j — 0, 1, 2) for scalar, vector, and tensor perturbations at given /. For n — 2 tensor 
perturbations do not exist. In this case the scalar and vector entries correspond to the QNMs of the Zerilli and 
Regge- Wheeler equations, which are known to be isospectral |21|. The isospectrality is broken for n > 2. This has 
been shown analytically by Kodama and Ishibashi |30j | and later verified numerically by Konoplya [33 | . The real and 
imaginary parts of scalar QNM frequencies at given n, I and j are smaller than those of vector QNMs, which are in 
turn smaller than those of tensor QNMs. Since scalar modes are the least damped, they are likely to dominate the 
gravitational radiation emission. 

As I grows, the isospectrality tends to be restored. In the eikonal limit / oo the centrifugal term of the potential 
dominates and is the same for scalar, vector, and tensor perturbations. In this limit, the QNM frequencies for all 
perturbations are 

The previous relation was derived in Rcf. 33] for multidimensional BH perturbations induced by a scalar field. (Notice 
that the normalization used in Ref. 33] is different from ours.) Here we have shown that it also holds for gravitational 
perturbations. Isospectrality of scalar and gravitational perturbations is a common feature of the eikonal limit and 
of the large-damping limit 23iE3 for any n. 



4.2. Multipolar components of the energy spectra 

The numerical integration described in Sec. Ogives the energy spectra of Figs, ^andl^l The spectra for n = 2 (top 
left panel in Fig. are in excellent agreement with those of Ref. [22]. The spectra are fiat in the region between the 
zero-frequency limit and a "cutoff" frequency ujc, beyond which they fall exponentially to zero. The cutoff frequency 
ujc is given by the fundamental QNM frequency to a good level of accuracy. This result can be understood in terms of 
gravitational-wave scattering from the potential barrier which surrounds the black hole, plays the role of the energy 
in the Schrodinger-like equation iQ. From Eq. (|15|) it follows that uj'^ — Vq at first order in the WKB approximation. 
Therefore, only the radiation with energy smaller than the peak of the potential is backscattered to infinity; radiation 
with larger frequency is exponentially suppressed. 

The gravitational emission of a two-particle hard collision in hi ghe r dimensions has been computed by Cardoso et al. 
[25I ] using techniques developed in four dimensions by Weinberg psj ] and later used by Smarr [23] . The main result of 
Ref. [2^ is that the spectra in n-l- 2 dimensions grow as w""^, thus the integrated spectra diverge as cj"^^. Physically 
meaningful results for the total energy can only be obtained by imposing some cutoff on the integrated spectra. Smarr 
[2^ ] first suggested to use the inverse horizon radius as a cutoff. The relativistic perturbative calculation in rt = 2 p3 ] 
shows that the cutoff frequency at fixed I is very close to the fundamental BH QNM. Therefore, the cutoff frequency 
should be given by some "weighted average" of the fundamental gravitational QNM frequencies E^l . 

Our results for the spectra and the QNMs confirm the above picture. Fig. ^ shows that all spectra go to zero as 
w — !■ 0. For Lu < LUc the spectrum at fixed / is 
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where fn,i is a constant that can be found by a fit of the spectra. For large I fn,i decays as 

/„. = fc„+2r3("+2)/4. (28) 

A fit of the numerical data gives — 2.25, fcg — 0.832, fcg = 0.184, and fcio = 0.040. Our result for n = 2 is consistent 
with that of Ref. [2^ . 

As conjectured in Ref. ^2^, all spectra have a maximum at some cutoff frequency uJc- This cutoff frequency is very 
close to the fundamental QNM frequency uiin for (scalar) gravitational perturbations with given I and n, which is 
marked by open circles in Fig. ^ and Fig. El The deviation between ujc and LOin is of order 10 % for low I, and 
decreases for large / (compare Fig. and Fig. El to the first column of Tables ITllV|l . The deviation is larger when the 
WKB method is least reliable, namely for I = 2 and n > 4. In these cases, the location of the peaks in the spectra 
can presumably be used as a more reliable estimate of the QNM frequency. The spectrum decays exponentially for 
oj > UJc with an n-dependent slope a„+2 (see Fig. jJl : 



duj 



(29) 



Thus the w-integrated multipolar contributions AEi at given I are finite. With our choice of units, Cardoso and 
Lemos find a4 = 13.5a (here a is a constant of order unity which cannot easily be determined because the spectra 
decay very quickly). Our numerical fits give ~ 15, in good agreement with their result. In higher dimensions the 
constants an+2 are comparatively easier to determine. Their values are ag — 5.5, as — 3.4, and aio — 2.3. It is not 
clear if there is any relation between this n-dependent slope and the late-time tail behavior predicted in Ref. . 

Fig- n shows that higher multipoles contribute more as n grows. This is evident when we look at the w-integrated 
multipolar components of the energy spectra of Fig.O The quadrupole {I = 2) is dominant only for n = 2 and n — A. 
For n = 6 and n — 8 the dominant multipoles are I = 4 and I = 6, respectively (see Table IV|) . This effect may be 
related to the appearance of a negative well in the scalar potentials for ^ = 2 and n > 4. It would be interesting to 
understand better the physical relation between the dominant multipole and the spacetime dimension. 

4.3. Total energy 

The total emitted energy is obtained by numerically integrating the results of the previous section over oj and 
summing the multipolar components. For large I the integrated energy in the multipole I can be fitted by 

AEi = an+2r''"+' , (30) 

where (og = 0.110,^6 = 1.69), (ag = 0.050,^8 = 1-64), and {aw = 0.022, 6io = 1.40) for n = 4, n = 6, and n = 8, 
respectively. The coefficients (a„+2, bn+2) have been obtained by fitting the data from ^ = 14 to Z = 20 and are weakly 
dependent on the chosen range of I. This variability affects our final results on the total energy within less than a few 
percent. 

Restoring the dependence on the BH horizon and on the conserved particle energy pq, the total emitted energy 

is 

^ ^ y ^ ^ - ^4"?^ , (31) 

Mbh 167r ^ Mbh 167r Mbh ^ ^ 

where E^l^^ is the "dimensionless" total energy (expressed in the units = 1 that we used in our numerical integra- 
tions). We obtained the integrated spectra numerically up to Z = 20 and extrapolated them for larger I using the fits 
in Eq. (|30|l . Results are presented in Table IVTI 

Following Ref. we estimate the gravitational energy loss for a collision of two particles with equal mass M by 
the replacement po ^ M, Mbh — ^ Aftot = 2M. For n — 2 this extrapolation gives results in good agreement with the 
perturbative shock-wave calculation of Ref. which considers two BHs of equal mass. An analogous extrapolation 
for n = 2 gives results in close agreement with the fully relativistic computation "iG] for a particle starting from infinity 
at rest. Therefore, we believe that our extrapolation should provide a qualitative but realistic estimate. The results 
for different values of n are given in the last column of Table IVT] The gravitational energy loss is 13%, ~ 10%, 7%, 
and 8% for n = 2 to n = 8, respectively. The result for n = 2 is in good agreement with previous estimates [13 (see 
the discussion in Ref. p^l. 
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4.4. Angular dependence 

The angular dependence of the radiation is obtained by evaluating numerically Eq. ljB17|l and Eq. (|B18|) . Fig. ^ 
shows the angular dependence of the total energy up to / = 15 for n = 2, n = 4 and n = 6. 

The angular distribution of the gravitational radiation in the BH frame goes to zero along the axis of the collision 
{6 = 0, tt) in any dimensions. Therefore, the gravitational emission is never back or forward scattered. In four 
dimensions the angular spectrum of the gravitational radiation increases rapidly at small 6 and becomes approximately 
flat at greater angles with a maximum in the direction orthogonal to the axis of the collision, before falling rapidly 
to zero for values of the angles close to tt. The angular distribution of the gravitational radiation for n > 2 is peaked 
at ±6 and tt ± 6*, where is a small angle. The difference between the behavior of the angular distribution in four 
dimensions and in higher dimensions has no evident physical reason. It would be interesting to further explore this 
point. 

5. CONCLUSION AND PERSPECTIVES 

In this paper we have computed the gravitational emission of a two-particle collision in an even £)-diniensional 
spacetime. We have presented the numerical results for _D = 4 to 10. The collision has been modelled as a massless 
test particle plunging into a BH with mass equal to the cm. energy of the event. 

According to our estimates, the total emitted energy in a head-on collision with particles of equal mass ranges from 
^ 13% (D = 4) to ~ 8% (D = 10). This shows that the loss in gravitational radiation is quite stable under variation 
of the spacetime dimension and slightly decreases for higher D. The result for D = A confirms previous numerical 
and analytical calculations |0| . 

Our result contrasts with the YN estimation for the initial mass of a BH in head-on collisions A possible 

explanation is that the junk emission is not wholly gravitational emission. The YN method predicts the mass within 
the apparent horizon to be ~ O.Tl-Ecm in four dimensions. If all the junk energy were gravitational radiation, this 
would amount to a total loss of around 30%. The disagreement is likely not due to numerical uncertainties or 
inaccurate approximations: the YN mass decreases for higher spacetime dimensions {'^ 0.71i?cm to ~ O.SSE'cm for 
D — A to D = 11), whereas the loss in gravitational radiation remains stable. Since both YN and our methods are 
purely gravitational, this "dark component" of the junk radiation should describe the by-products of the collision. 
According to this picture, ~ 60% of the cm. energy in ten dimensions is trapped inside the horizon, ~ 10% is emitted 
in gravitational radiation, and ^ 30% goes into particle by-products in the final state. These could be the carriers of 
the initial charge in a collision between charged particles. For _D > 4 a fraction of the by-products may be emitted 
into the bulk. 

Let us conclude by briefly discussing the phenomenological consequences of these results for BH formation at the 
TeV scale. Although uncertainties may affect the numerical estimates, different approaches now confirm that some of 
the initial cm. energy is not trapped inside the BH horizon. For head-on collisions in _D = 10, for example, this junk 
energy ranges from ^ 10% (optimistic value ~ our result) to ^ 40% (pessimistic value - YN result). Hence, the initial 
mass of the BH formed in the collision could be considerably smaller than the cm. energy. The experimental signatures 
of BH production at particle colliders and in ultrahigh energy cosmic ray events strongly depend on the initial BH 
mass. The total multiplicity of the Hawking phase in ten dimensions could be almost halved in the pessimistic case, 
leading to a greater average energy of the emitted quanta. 

A thorough investigation of the effects of energy loss in TeV-scale BH production is undoubtedly worth pursu- 
ing. Future research should focus on the extension of the above results to spacetimes with odd dimensions and to 
gravitational events with different geometries. BHs produced in colliders, for instance, possess nonvanishing angular 
momentum. Rotating BHs are expected to lose more energy in gravitational waves than Schwarzschild BHs of equal 
mass. A larger gravitational emission is also expected for non-spherically symmetric BHs. This is particularly relevant 
when the compactified space is asymmetric, and some of the extra dimensions have size of order of the fundamental 
gravitational scale. It would be extremely important to quantify these differences. 
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APPENDIX A: THE SOURCE TERM 



In this Appendix we derive the source term of the KI equation that describes the radial plunge of a massless particle 
into the {n + 2)-diinensional BH. The perturbation of the stress-energy tensor is 



Po 



'dX'dX 



(Al) 



where po is the conserved energy of the particle. The only non-vanishing components of the particle velocity are u* 
and . Thus the source excites only scalar perturbations. Following the notations of Ref. ,22i|, Eq. (|A1() reads 



5Tny = 



TabS 












(A2) 



where S are the scalar harmonics and Tab are the non-vanishing gauge-invariant perturbations of the stress-energy 
tensor. The BH+source system is symmetric under rotation of the {n — l)-sphere S'^~^ Consequently, the 

harmonic decomposition of the fields contains only harmonics invariant under S"""^. We can write the metric of S*" 
as 



+sin2 OdVln-Ml, ■ 



(A3) 



and choose the trajectory of the test particle to be = 0. The harmonics which are invariant under S" ^ do not 
depend on . . . , (/)„. The scalar harmonics on 5" belong to the representations D^'^ of SO{n + 1) 



(A4) 



Each harmonic is labelled by the index I denoting its representation and by additional indices in the representation. 
We fix a particular element of each representation D^'^ (the singlet under SO{n — 1)) by requiring the harmonics to 
be invariant under 5""^. Therefore, the harmonics in the expansion of the perturbations depend only on / and on 
the dimension n of the sphere. A {n + 2)-dimensional scalar field can then be expanded as [ig 

ct>{t, r, 0, </>!,..., 0„_i ) = ^ 0; (i, r)5("') {9) , 



where 5^"') {0) satisfy 



and 



The solution of Eq. (|A6|) is 



(A5) 

(A6) 
(A7) 

(A8) 



where ^^^^{9) are Gegenbauer polynomials ^3] and iCf"') are normalization factors. Using Eq. (jA7|) we have 

-1/2 

(A9) 



T{l + n- 1) 



+ f-5)r(f-i) r(; + i) 



The scalar harmonics for the source are obtained setting 9 = 0: 

^ ^ (n-2)!/! 



For a massive particle in radial geodesic motion 



dt 
dr 



fir) 



(AlO) 



(All) 
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From equations (|A11|I and l@J it follows that 

r,{r)^-t{r), (A12) 
where we have set to zero the integration constant. The stress-energy tensor perturbations are 

<5r„(i,r,0,0i,...) = -po.rV-",5(t-t(r)),5"(r!„), 

5Ttr{t,r,6,(^^,...) - -porV-"<5(t-t(r)),5"(r!„). (A13) 

{5Ttt does not contribute to the source of the KI equation.) Integrating Eqs. (|A13|I on S'" and applying a Fourier- 
transform, the gauge-invariant perturbations are 

piujt(r) 

Trr{cO,r) = -po/-^-" — ^5^"') (0 = 0) , 

V ZTT 

i(^t(r) 

Ttr{io,r) = -po/-V-"^^5("')(0 = O). (A14) 

The source term S*^"^ for scalar gravitational perturbations is obtained in terms of Tab by substituting Eqs. HA14|I in 
Eq. (5.44) of Ref. H^, where 



Sab - Snr^'-^Tab, Sa = St^A^ Ja^O. (A15) 



The result is: 



^(.) ^ ^-....8V4rr2(/-i)(/ + 2)(.-i) 

i^r [(; + 2)(;-l)r + 3]' 

c,(4) ^ .-i.., 16V(2/ + 3)A2 [(^^ + 30(r^ - r^) + 5 - 4r^ - r^] 

\/2^ic^r3 [a + 4)a-l)r3 + 10]2 ' ^ ' 



Si6) ^ ^-i.r, 6V(2^ + 5)A4 [jP + 50(ri" - r^) + 14 - Gr^" - Sr^] ^ ^^^^^ 



[(; + 6)(;-l)r5 + 21] 
^(8) ^ -i.., Wi2l + 7)A6 [(P + 70(ri4 - r^) + 27 - Sr^^ - 19r^] 

V3^iu;r5 [(Z-h8)G-l)r7 + 36]2 ' ^ 

where Aj, = [{I + k){l + k — 1) . . . {I + 1)]. The BH horizon r/^ and the particle energy pq have been set equal to one. 
Equation l|IT6|) coincides with Eq. (8) of Ref. 

APPENDIX B: ENERGY AND ANGULAR DISTRIBUTION 

We derive here the formulae for the energy and the angular dependence of t he g ravitational emission. Gravitational 
waves in a (n -I- 2)-dimensional spacetime behave asymptotically as ~ r~"/^ [23 and possess n{n + l)/2 — 1 degrees 
of freedom. The transverse-traceless (TT) gauge is defined by Sgab = 0, Sgai = 0, and SgijY'' = 0. These conditions 
can be chosen by imposing the harmonic gauge in the wave zone and using the remaining gauge freedom to constrain 
Sgtt = 0, Sgu = 0, and J'■^Sg^j =0. 

We separate the angular part of the perturbations using tensor spherical harmonics, following Ref. [30|- In the TT- 
gauge, the only nonvanishing term in the decomposition is the Ht component. Hence, the gauge invariant quantities 
for scalar perturbations {F and Fab) depend only on Hx - 

F=-HT + -{D''r)Xa, Fab = DaXb + DbXa, (Bl) 

n r 

where 

Xa = ^DaHr . (B2) 
The scalar perturbation <I> is written in terms of the gauge invariant quantities as 

^ nZ-r{X + Y) 

W2-1[k2 + l).x/2] ' ^ ' 
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where 



Setting 



the asymptotic behavior of <f> is 



X + Y = -2nr''-'^F , Z = —r'^~'^F% . (B4) 

iw 



Ht - -^e^-' , (B5) 



OA 

Hm $ = ^e'""'* . (B6) 



The asymptotic behavior of hij in the TT-gauge is 

= ^-"HtS.^ - ^^^^ - ^ {d^D.S + ^7..5) . (B7) 
The energy-niomentum pseudotensor does not depend on the spacetime dimension and is given by pH 

'^^ =<t™>-:^</i^^/i^^'^">, (B8) 



dSdt 32n 

where hjj^ are metric perturbations in the time domain. Using Parscval's theorem, the energy-momentum pseudoten- 
sor in the frequency domain is 

'^^ <Wh*^^^'>, (B9) 



where hjj'" are now the metric perturbations in the frequency domain. Substituting Eq. ljB7(l in Eq. ljB9(l we get 

Integrating on the sphere dS — r"d^ln and using the relations 

D.D'S = ~K,^S , [Dj,D,]V^ = RuV' = [n~ l)V, , (Bll) 

where is a generic vector (here and in the foUowing, indices are raised and lowered with 7*-' ), we find the "two-sided" 
power spectrum 

aUJ two-sided J dbdui dZTT Tl 

The "one-sided" spectrum in Eq. (|21|l is obtained by multiplying Eq. ljB12p by two. In the four-dimensional limit the 
one-sided spectrum is 



dE 

duj 



Equation (|B13p is Zerilli's formula for the l-th multipole component of the energy spectrum in four-dimensions. The 
total energy spectrum is given by the sum over the multipoles. 
For each I the energy spectrum is 

Substituting Eq. ljB12p in the previous equation, we find 

dEi _dEi n 1 
drtdw dw n — 1 k'^{k'^ — n) 



D.DjS + — 7y5 D'D^S* + —7*^5* . (B15) 
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The angular dependence for the Z-th multipole is obtained by integrating over frequency '2^ . The result is 
dEi . _ n 1 

where 



AK^) = ^rT — T f^'5.«« + -^'5) ■ (B17) 

For n = 2 Eq. ljB17p reduces to the known result in four dimensions The angular dependence is obtained by 

summing over the multipoles: 

The result truncated to Imax = 15 is shown in the left panel of Fig. 01 The curve corresponding to ?i = 2 shows good 
agreement with Fig. 3 of Ref. . 
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CO 



FIG. 1: Multipolar components of the energy spectra up to Z = 10 for n = 2, n = 4 (top left and top right panels), n = 6 
and n = 8 (bottom left and bottom right panels) in units th = 1- Open circles mark the real part of the fundamental scalar 
gravitational quasinormal frequency, LOin, for the given I and n. 




FIG. 2: Multipolar components of the energy spectra at fixed I for different values of n in units Vh = 1. The left panel 
corresponds to Z = 2 and the right panel corresponds to Z = 3. Open circles mark the real part of the fundamental scalar 
gravitational quasinormal frequency, uJin, for the given Z and n. 
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FIG. 3; The integrated energy AEi as a function of I for different values of n. The dominant muhipolar component is / = 2 
only for n < 6; this is probably related to the appearance of a negative well in the scalar potentials for / — 2 and n > 4. The 
dominant multipole is Z = 4 (6) for n = 6 (8) (see Table 0. 




FIG. 4: Angular dependence of the radiation for n = 2, 4 and 6, summing all multipoles up to / = 15 in units rh ~ 1. The 
angular distribution for n = 8 is not shown. The latter is even more peaked in a narrow region around ^ = since also the 
multipoles with / > 15 contribute significantly to the radiation. 
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TABLE I: QNMs for n = 2. The first three quasinormal frequencies for scalar and vector perturbations are listed from left to 
right. The scalar modes and the vector modes are isospectral in four dimensions. 
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TABLE 11: QNMs for n = 4. The first three quasinormal frequencies for scalar, vector, and tensor perturbations are listed 
from left to right. 
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6.852-1 


4591 


6.725-2.4441 


6.992-0.4891 


6.928-1 


4701 


6.801-2 


4631 


7.138-0.4951 


7.074-1 


4881 


6.948-2.4921 


12 


7 


504.0.4871 


7.444-1 


4631 


7.326-2.4481 


7.574-0.4901 


7.514-1 


4721 


7.396-2 


4641 


7.708-0.4951 


7.649-1 


4871 


7.532-2.4891 


13 


8 


088-0.4881 


8.033-1 


4651 


7.923-2.4521 


8.153-0.4901 


8.098-1 


4731 


7.988-2 


4651 


8.279-0.4951 


8.224-1 


4861 


8.115-2.4871 


14 


8 


671-0.4891 


8.619-1 


4681 


8.517-2.4541 


8.732-0.4911 


8.680-1 


4741 


8.577-2 


4651 


8.849-0.4951 


8.798-1 


4861 


8.696-2.4851 


15 


9 


253-0.4891 


9.204-1 


4691 


9.108-2.4561 


9.309-0.4911 


9.261-1 


4751 


9.164-2 


4661 


9.420-0.4951 


9.371-1 


4851 


9.275-2.4831 
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TABLE III: QNMs for n = 6. The first three quasinormal frequencies for scalar, vector, and tensor perturbations are listed 
from left to right. The numbers in italic indicate that the potential at the given I is not everywhere positive definite. The square 
brackets indicate that the potential has two scattering peaks. 



n=4 


Scalar modes 


Vector modes 


Tensor modes 


1 


j = 





j = 


1 


3 = 


2 


3 - 







3 = 


1 


3 - 




2 


j = 


j = 1 


3 - 




2 


2 


[1. 778-0. 57n] 


[1.289-1 


77 (A] 


l0.395-3.201i] 


2.388 


-0.72(K 


1.831-2 


237\ 


0.825 


-4 


OOli 


2.956-0.7511 


2.365-2.3571 


1.339- 


4 


245i 


3 


2.604-0 


628i 


2.198-1 


9161 


1.403-3 


3551 


3.102- 





7151 


2.660-2 


1911 


1.814- 


3 


8331 


3.623-0.7471 


3.181-2.2941 


2.351- 


4 


0121 


4 


3.401-0 


645i 


3.050-1 


9581 


2.346-3 


3751 


3.815- 





7121 


3.450-2 


1651 


2.730- 


3 


7311 


4.282-0.7441 


3.926-2.2641 


3.235- 


3 


8951 


5 


4.174-0 


660i 


3.875-1 


9971 


3.270-3 


4031 


4.522- 





7121 


4.213-2 


1561 


3.595- 


3 


6781 


4.940-0.7411 


4.640-2.2471 


4.047- 


3 


8301 


6 


4.923-0 


675i 


4.665-2 


0371 


4.144-3 


4491 


5.222- 





7141 


4.954-2 


156i 


4.418- 


3 


6541 


5.598-0.740i 


5.337-2.236i 


4.818- 


3 


789i 


7 


5.653-0 


687i 


5.425-2 


0701 


4.967-3 


4921 


5.915- 





7161 


5.679-2 


1601 


5.207- 


3 


6451 


6.255-0.7391 


6.024-2.2291 


5.563- 


3 


7631 


8 


6.369-0 


695i 


6.164-2 


0951 


5.753-3 


525i 


6.602- 





7191 


6.392-2 


1641 


5.969- 


3 


6431 


6.913-0.7381 


6.705-2.2241 


6.290- 


3 


745i 


9 


7.075-0 


702i 


6.888-2 


1131 


6.515-3 


5501 


7.285- 





7211 


7.094-2 


1691 


6.712- 


3 


6441 


7.570-0.7381 


7.382-2.2211 


7.004- 


3 


7321 


10 


7.772-0 


7071 


7.602-2 


1281 


7.259-3 


5701 


7.964- 





7221 


7.790-2 


1731 


7.441- 


3 


6461 


8.228-0.7371 


8.055-2.2181 


7.709- 


3 


7221 


11 


8.464-0 


7111 


8.306-2 


1391 


7.989-3 


5851 


8.640- 





7241 


8.480-2 


1771 


8.158- 


3 


6481 


8.885-0.7371 


8.726-2.2161 


8.406- 


3 


7151 


12 


9.151-0 


715i 


9.004-2 


1481 


8.709-3 


5981 


9.314- 





7251 


9.165-2 


1801 


8.867- 


3 


6501 


9.543-0.7371 


9.395-2.2151 


9.098- 


3 


7091 


13 


9.834-0 


7171 


9.697-2 


1561 


9.421-3 


6071 


9.986- 





7261 


9.847-2 


1831 


9.569- 


3 


653i 


10.200-0.7371 


10.062-2.2141 


9.785- 


3 


705i 


14 


10.51-0 


720i 


10.39-2 


1621 


10.13-3 


6161 


10.66- 





7271 


10.53-2 


1851 


10.27- 


3 


6551 


10.86-0.7371 


10. 73-2. 2131 


10.47- 


3 


7011 


15 


11.19-0 


7211 


11.07-2 


1671 


10.83-3 


6221 


11.32- 





7281 


11.20-2 


1871 


10.96- 


3 


6571 


11.52-0.7361 


11.39-2.2121 


11.15- 


3 


6981 



TABLE IV: QNMs for n = 8. The first three quasinormal frequencies for scalar, vector, and tensor perturbations are listed 
from left to right. The rmmbcrs in Italic Indicate that the potential at the given I is not everywhere positive definite. The square 
brackets indicate that the potential has two scattering peaks. 



n=4 


Scalar modes 


Vector modes 


Tensor modes 


1 


3 = 


i = 1 


3 = 2 


3 = 


i = i 


3 = 2 


.1 = 


.i = i 


.1 = 2 


1 = 2 


[2.513-0.744ij 


[1. 686-2. 299iJ 


[0.159-4. 345ij 


3.261-0.924\ 


2. 335-2. 85 li 


0.598-5.2871 


3.886-0.9591 


2.765-2.9881 


0.706-5.720i 


1 = 3 


[3.388-0.8121] 


[2.696-2.4611] 


[1.277-4.431i] 


4.017-0.9231 


3.269-2.804\ 


1.747-5.0161 


4.618-0.959i 


3.806-2.917i 


2.141-5.2411 


1 = 4 


4.223-0.8411 


3.631-2.5321 


2.367-4.4201 


4.775-0.9201 


4.147-2.7771 


2.824-4.8401 


5.336-0.9551 


4.691-2.8851 


3.331-5.0181 


1 = 5 


5.042-0.8551 


4.524-2.5681 


3.407-4.4011 


5.531-0.9181 


4.991-2.7621 


3.840-4.7341 


6.049-0.9511 


5.507-2.8661 


4.360-4.9041 


1 = 6 


5.848-0.865i 


5.390-2.595i 


4.403-4.399i 


6.283-0.917i 


5.810-2.7571 


4.802-4.676i 


6.761-0.9491 


6.291-2.854i 


5.297-4.838i 


1 = 7 


6.640-0.8741 


6.231-2.6221 


5.357-4.4151 


7.030-0.9181 


6.610-2.7571 


5.719-4.6461 


7.473-0.9471 


7.056-2.8461 


6.178-4.7981 


1 = 8 


7.420-0.8831 


7.052-2.6471 


6.270-4.4411 


7.774-0.9201 


7.396-2.7591 


6.599-4.6341 


8.184-0.9461 


7.808-2.8411 


7.022-4.7721 


1 = 9 


8.191-0.8901 


7.855-2.6691 


7.149-4.4691 


8.514-0.9211 


8.167-2.7641 


7.450-4.6301 


8.895-0.9451 


8.553-2.8371 


7.841-4.7551 


1 = 10 


8.953-0.896i 


8.645-2.688i 


7.999-4.495i 


9.250-0.923i 


8.935-2.7681 


8.278-4.630i 


9.606-0.944i 


9.292-2.834i 


8.640-4.743i 


1 = 11 


9.709-0.9021 


9.423-2.7051 


8.829-4.5191 


9.984-0.9241 


9.692-2.7731 


9.088-4.6331 


10.32-0.944i 


10.03-2. 8321 


9.426-4.7351 


1 = 12 


10.46-0.9061 


10.19-2.7181 


9.641-4.5391 


10.71-0.9261 


10.44-2.7771 


9.884-4.6371 


11.03-0.9431 


10.76-2.8311 


10.20-4.7281 


1 = 13 


11.20-0.910i 


10.95-2.730i 


10.44-4.556i 


11.44-0.927i 


11.19-2.7811 


10.67-4.642i 


11.74-0.943i 


11.49-2.829i 


10.97-4.724i 


1 = 14 


11.95-0.9141 


11.71-2.7401 


11.23-4.5711 


12.17-0.9281 


11.93-2.7851 


11.44-4.6461 


12.45-0.9431 


12.21-2.8281 


11.73-4.7201 


1 = 15 


12.68-0.9161 


12.46-2.7481 


12.01-4.5841 


12.90-0.9291 


12.67-2.7881 


12.21-4.6501 


13.16-0.9431 


12.94-2.8281 


12.48-4.7181 
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TABLE V: Multipolar contributions to the total energy for different n in units = 1. 



I 


n — 2 


n = 4 


n = 6 


n = 8 


2 


0.1845 


O.lSQo- 


1 





194o- 


2 


0.269c- 


3 


3 


0.0855 


0.120c- 


1 





238o- 


2 


0.653O- 


3 


4 


0.0500 


0.086e- 


1 





241e- 


2 


0.983o- 


3 


5 


0.0329 


0.064e- 


1 





224o- 


2 


1.187o- 


3 


6 


0.0234 


0.050o- 


1 





199o- 


2 


1.258c- 


3 


7 


0.0175 


0.039e- 


1 





172e- 


2 


1.225e- 


3 


8 


0.0136 


0.032o- 


1 





149o- 


2 


I.I3O0- 


3 


9 


0.0109 


0.027O- 


1 





128o- 


2 


1.009c- 


3 


10 


0.0089 


0.022O- 


1 





lllo- 


2 


0.8880- 


3 


11 


0.0074 


0.019c- 


1 





096o- 


2 


0.780c- 


3 


12 


0.0063 


O.OI60- 


1 





084o- 


2 


0.688c- 


3 


13 


0.0054 


0.014O- 


1 





075e- 


2 


0.612c- 


3 


14 


0.0047 


0.013e- 


1 





O660- 


2 


0.549c- 


3 


15 


0.0041 


O.Olle- 


1 





059e- 


2 


0.497e- 


3 



TABLE VI: Total energy for different spacetime dimensions. From left to right, the columns give the spacetime dimension 
D = n + 2, the factor n^„/167r, the total energy E^l^^ in units ru = 1, the rescaled total energy £^^^ = nAnE^I^^ /IGir and the 
gravitational energy loss (see text). 



D 






''tot 


Energy loss 


4 


1/2 


0.52 


0.26 


13% 


6 


27r/3 


0.095 


0.20 


10% 


8 


27rV5 


0.034 


0.13 


7% 


10 


167r^/105 


0.032 


0.15 


8% 



